SIMILAR SUBLATTICES AND COINCIDENCE ROTATIONS 
OF THE ROOT LATTICE A 4 AND ITS DUAL 
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Abstract. A natural way to describe the Penrose tiling employs the projection method 
on the basis of the root lattice A4 or its dual. Properties of these lattices are thus related 
to properties of the Penrose tiling. Moreover, the root lattice A4 appears in various other 
contexts such as sphere packings, efficient coding schemes and lattice quantizers. 

Here, the lattice A4 is considered within the icosian ring, whose rich arithmetic structure 
leads to parametrisations of the similar sublattices and the coincidence rotations of A4 and 
its dual lattice. These parametrisations, both in terms of a single icosian, imply an index 
formula for the corresponding sublattices. The results are encapsulated in Dirichlet series 
generating functions. For every index, they provide the number of distinct similar sublattices 
as well as the number of coincidence rotations of A4 and its dual. 



1. Introduction and generalities 

A lattice r in R d , denoted by r = (61, ... , bd)z, consists of all integer linear combinations 
of a basis {61, . . . , b^} of R d . In crystallography, coincidence site lattices of a lattice r are 
used to describe and understand grain boundaries; compare [2] and references given there. 
Mathematically, a coincidence site lattice (CSL) of a lattice r is defined as r n Rr, where R 
is an isometry and the corresponding coincidence index 

E(R) = [r ■. (rn Rr)] 

is finite, i.e. R and Rr share a common sublattice. Denote the orthogonal group by 0(d,R) 
and define the set of all coincidence isometries as 

OC(r) := {R G 0(d,R) I E{R) < 00}. 

According to [2], it is a subgroup of 0(d, R). OC(-T) and SOC(r'), which is the subgroup of 
all coincidence rotations, can be characterised as subgroups of another important group of 
isometries, namely 

OS(r) := {R G 0(d, R) I aRr C T, a G R, a > 0}. 

This group consists of all isometries that emerge form similarity mappings of the lattice r into 
itself; see [3] for details. Its subgroup of rotations, denoted by SOS(r), contains SOC(r) as a 
normal subgroup; see [10] for a detailed analysis of their relation. A sublattice of r of the form 
aRr is called a similar sublattice (SSL) of r and has obviously the index [r : aRr] = a d . If 
the lattice r possesses a rich point symmetry, there are usually interesting SSLs, beyond the 
trivial ones of the form mr with m G N. For R G OS(.T) define the denominator relative to 
the lattice r as 

den r (i?) := min{a G R | a > 0, aRr C T} . 
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As R is an isometry, one always has den r (R) > 1, and from den r (R)Rr C r one concludes 

that (den r (i?)) is an integer. Consequently, den r (R) is either a positive integer or an 
irrational number, but still algebraic and of degree < d. Moreover, one has 

{a G R | a > 0, aRr Cf} = den r (R) N . 

This leads to the following characterisation of OC(/ 1 ) within OS(-T), compare [3] for details: 

(1) OC(r) = {i? G OS(r) | den r (#) G N}. 

When a lattice J 1 C M rf is given, one is interested in the number of distinct SSLs of r of index 
n as well as in the number of distinct CSLs of r of index n. If these arithmetic functions are 
multiplicative, they are often encapsulated into Dirichlet series generating functions, because 
of their Euler product decomposition. A detailed introduction to arithmetic functions, the 
corresponding Dirichlet series and Euler products can be found in pQ. For many lattices in 
d < 4 the arithmetic functions which count the number of distinct SSLs and CSLs of each 
index have already been derived; see for instance [21 El [7]. One lattice for which this problem 
has not been solved completely is the root lattice A4. This lattice is of particular interest, 
because it forms the natural setting, in the sense of a minimal embedding, for the description 
of the Penrose tiling as a cut and project set, see for example [5j. Properties of the root lattice 
A4 are thus directly related to the Penrose tiling. Some of the various other applications of 
the root lattice A4 are described in [9]. 

In this paper recent results about the similar sublattices and coincidence rotations of the 
root lattice A4 are reviewed and extended to its dual lattice A\; details can be found in 
[U [3]. In the following sections we first describe the root lattice A4 and its dual A\ in a 
suitable setting. Then we derive their Dirichlet series generating functions for the number of 
distinct SSLs of index n, which turn out to be the same for A4 and A\. For these lattices it 
is particularly difficult to determine the Dirichlet series generation functions for the number 
of distinct CSLs of index n. Therefore, we consider the slightly simpler problem to derive the 
generating functions for the number of distinct coincidence rotation of index n, which turn 
out to coincide, too. 

2. The root lattice A$ 

The root lattice A4 is usually defined as 

A A := {(xi,...,x 5 ) G Z 5 I xi + ... + s 5 = 0}, 

which lies in a 4— dimensional hyperplane of IR 5 . One description of A4 as a lattice in M 4 is 
given by 

L := ((1,0,0,0), i(-l, 1,1,1), (0,-1, 0,0), i(0,l,r-l,-r)) z 

where r = (l + \/5)/2 is the golden ratio. Relative to the inner product ti(xy) = 2(x\y), where 
(x\y) denotes the standard Euclidean inner product, this lattice is the root lattice A^; see 
[HI SI El for details. This particular description of the root lattice A± in R 4 is very convenient 
for our problem, as it enables us to use the arithmetic of the quaternion algebra H(Q(\/5 )); 
see [12] for a detailed introduction to Hamilton's quaternions. For brevity we use from now 
on the notation 

K := Q(V5 ) = {q + rVE | q, r G Q}, 
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which is a quadratic number field. The algebra M(K) is explicitly given as 

M(K) = K®iK®]K®kK, 

where the generating elements satisfy Hamilton's relations i 2 = j 2 = k 2 = ijk = —1. It is 
equipped with a conjugation 7 which is the unique mapping that fixes the elements of the 
centre of the algebra K and reverses the sign on its complement. If we write 

q = (a,b,c,d) = a + ib + jc + kd, this means q = (a,—b, — c, —d) . 

The reduced norm and trace in M(K) are defined by 

m(q) := qq = \q\ 2 and tr(g) := q + q, 

where we canonically identify an element a G K with the quaternion (a, 0,0, 0). For any 
q G M(K), \q\ is its Euclidean length, which need not be an element of K. Nevertheless, 
one has \rs\ = \r\\s\ for arbitrary r, s G M(K). Due to the geometric meaning, we use the 
notations \q\ 2 and nr(g) in parallel. An element q G M(K) is called integral when both m(q) 
and tr(q) are elements of 

:= Z[r] := {m + nr | m,n 6 Z}, 

which is the ring of integers of the quadratic field K. 

The icosian ring I consists of all linear combinations of the vectors 

(l,0,0,0),(0,l,0,0),i(l,l,l,l),i(l-r,r,0,l) 

with coefficients in o. The elements of I are called icosians. It is a remarkable object; see for 
example [151 El El Sj and references given there. In particular, I is a maximal order in M(K), 
all elements of I are integral in M(K) and the lattice L is contained in I. Since has class 

number 1, compare |15} I16j. all ideals of I are principal. The detailed arithmetic structure of 
I is the key to the characterisation of the similar sublattices and coincidence rotations [U [3] 
for L. What is more, one significantly profits from another map, called the twist map in [4], 
which is an involution of the second kind for M(K); see pj] for details. If q = (a, b, c, d), it is 
defined by the mapping 

q^q:= (a' ,b' , d' , c) , 

where ' denotes the algebraic conjugation in K, as defined by the mapping \/5 i— > — The 
algebraic conjugation in K is also needed to define the absolute norm on K, via N(a) = \aa'\ . 
For the various properties of the twist map, we refer the reader to [4] and references therein. 
The most important ones in our present context are summarised in the following Lemma, 
whose proof can be found in [3]. It describes the relations between L and I. 

Lemma 2.1. Within M(K), one has 1 = 1 and 

L = {x £ I | x = x} = {x + x | x G 1} = </> + (I), 
where the Q-linear mapping : M(K) — > M(K), is defined by 4> + (x) = x + x. 

The dual A\ of the the root lattice A4, here in form of the dual of the lattice L, is given by 

L* := {x G M 4 I {x\y) G Z for all y G L}. 
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3. Similar Sublattices 

We are interested in the SSLs of the root lattice A4 and its dual lattice A\. According to 
[1], there is an index preserving bijection between the SSLs of A\ and A4, as well as between 
the SSLs of A4 and L. Therefore, it is sufficient to concentrate on the SSLs of the lattice L. 
For convenience, we define M(K)' = M(K) \ {0}. Proofs for this section can be found in [1]. 

Lemma 3.1. All SSLs of the lattice L are images of L under orientation preserving mappings 
of the form x \—* pxq, with p,q £ M(K)' . 

This characterisation tells us that we only need to select an appropriate subset of 

M(K)' x M(K)' in order to reach all SSLs of L. A first step is provided by the following 

observation. 

Lemma 3.2. If p £ I, pLp is an SSL of L. 

On the other hand there is the following strengthening of Lemma 13.11 

Proposition 3.3. If pLq C L withp,q £ M(K)', there is an a £ Q such that q = ap . 

To improve this characterisation we need the following definitions. An element p £ I is 
called I-primitive when ap £ I, with a £ K, is only possible with a £ 0. Similarly, a sublattice 
A of L is called L -primitive when a A C L, with a £ Q, implies a £ Z. Whenever the context 
is clear, we simply use the term "primitive" in both cases. 

Corollary 3.4. All SSLs of the lattice L are images of mappings of the form x 1— ► apxp with 
p £ I primitive and a £ Q. 

We now need to understand how an SSL of L of the form pLp with an I-primitive quaternion 
relates to the primitive sublattices of L. 

Proposition 3.5. If p £ I is I-primitive, pLp is an L-primitive sublattice of L. 

Combining the results of Corollary 13.41 and Proposition 13.51 we obtain the following impor- 
tant observation. 

Proposition 3.6. A similar sublattice of L is primitive if and only if it is of the form pLp 
with a primitive element p £ I. Moreover, all SSLs of L are of the form qLq with q £l. 

The next step is to find a suitable bijection that permits us to count the primitive SSLs of 
L of a given index. Recall from |14t [8] that the unit group of I has the form 

I x = {x £ I I N(nr(a?)) = ±1}, 

which leads to the following equivalence. 

Lemma 3.7. For p £ I, one has pLp = L if and only if p £ I x . 

Observe now that pi = ql with p, q £ I holds if and only if q _1 p £ I x . The relevance of 
this fact in our context comes from the observation that 

[I: pi] = N(H 4 ) = [L : pLp] , 

where the index of pi in I follows from the determinant formula in [H Fact 3] followed by 
taking the norm in Z[r]. This means that pi < — ► pLp describes an index preserving bijection 
between primitive right ideals of I (meaning right ideals ql with q £l primitive) and primitive 
SSLs of L. Due to the definition of L-primitivity, a general SSL can be described as an integer 
multiple of a primitive SSL. This leads together with the observation that all possible indices 
of SSLs are squares, to the following central result. 
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Theorem 3.8. The number of SSLs of a given index m is the same for the lattices A\,A^ 
and L. The possible indices are the squares of non-zero integers of the form k 2 + k£ — I 2 . If 
fssh( m ) denotes the number of SSLs of index m 2 , the corresponding Dirichlet series generating 
function reads 



Dssl(s) := 



/ssl(to) 



m=l 



m 



2s 



C(4s) 



CiOO 



where ((s) = *s the Riemann zeta function, 

p 

cm=j^ n 



i 



l-5~ s 11 1-p" 
p=±2(5) 



27 11 (i-p- s y 2 

p=±l(5) 



is the Dedekind zeta function of the quadratic field K , and Cj(s) = Cjr(2s) Ck(^ s ~ -0 denotes 
the Dirichlet series for the right ideals of I, which is the zeta function of I. 

Inserting the Euler products of £(s) and Cft-( s )> one finds the expansion of the Dirichlet 
series -Dssl(s) as an Euler product 



Dssl(s) 



(l-5- 2s )(l-5 1 - 2s ) 



n 



n i 



p=±l (5) p=±2 (5) 

Consequently, the arithmetic function fssh( m ) is multiplicative and therefore completely 
specified by its values at prime powers p r with r > 1. These are given by 

( 5 r+1 -l 



fssh(p r 



2 (l-p r+1 )-(r+l)(l-p 2 )p' 

(1-P) 2 
2-p r -p r + 2 
T^p 1 ' 

0, 



p = 5, 
p = ±l(5), 
p = ±2 (5), r even, 
p = ±2 (5), r odd. 



The first few terms of the Dirichlet series thus read 



n ( e\ 1_|_6_i_6ill_i_ 24 i 26 , 40 , 36 , 31 , 60 , 64 , 66 , 

-L^SSLyS ) r -\- ^ s -|- ~r g2s ~r \\2s ~r \q'2s ~t \g'2s ~r 20' 2s ' 25 2s 29 2s 312s ~r 3g2s "T ■ • ■ 

where the denominators are the squares of the integers previously identified in [9] . 

4. Coincidences and Rotations 

According to [21 Th. 2.2], 00(^4) = OC(A|) and the coincidence index of any conincidence 
isometry is the same for both lattices. Therefore, it is sufficient to investigate the CSLs of 
the lattice A±. The best representation of this lattice is again the lattice L. First of all the 
investigation can be restricted to coincidence rotations without missing any CSLs, because 
L = L, i.e. any orientation reversing operation can be obtained from an orientation preserving 
one after applying conjugation first. According to Proposition 13. 6\ a given SSL of L is of the 
form qLq, with q £ I. The corresponding rotation is given by the mapping x ^ yij- qxq. 
Of course, many different icosians q result in the same rotation. Our aim is to restrict q 
to suitable subsets of I without missing any rotation. In general, SOC(L) and SOS(L) are 
related according to Eq. (fTJ). Therefore, we have to identify the elements of SOC(L) within 
SOS(L). We would like to refer the reader to [3] for the proofs of this section. 
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Lemma 4.1. Let 7^ q G I be an arbitrary icosian. The lattice L n tAt gLg is a CSX o/L i/ 
and onZy i/ |g<f| G N. 

Let us call an icosian q G I admissible when |gg| G N. As nr(g) = nr(g)', the admissibility 
of q implies that N(nr(g)) is a square in N. With this definition the CSLs of L can be 
characterised as follows. 

Theorem 4.2. The CSLs of L are precisely the lattices of the form L n qLq with q G I 
primitive and admissible. 

This is the first step to define a bijection between certain primitive right ideals ql of the 
icosian ring and the CSLs of L. The next step in this direction is the following Lemma. 

Lemma 4.3. Let r, s G I &e primitive and admissible icosians, with rl = si. Then, one has 

rLf sLs 

Ln - — r = Ln- — -. 

\rf\ \ss\ 

For our further discussion we need to replace the primitive and admissible icosian p, and 
with it p, by certain o-multiples, such that their norms have the same prime divisors in 0. 
In view of the form of the rotation x 1— > j^pxp, this is actually rather natural because it 
restores some kind of symmetry of the expressions in relation to the two quaternions involved. 
For a primitive and admissible icosian p £l we choose explicitly 



/ lcm(nr(g),nr(g)) 

( ) aq = \ Mti G ' 

where we assume a suitable standardisation for the 1cm of two elements of 0; see [3] for details. 
Moreover, we have the relation a- = oT q = ct' q . The icosian a q q is called the extension of the 
primitive admissible element q G I, and (a q q, a' q q) the corresponding extension pair. Since 
a q and a' q are central, the extension does not change the rotation, i.e. 

qxq _ q a xq a 

\qq\ \Qaq a \ 

holds for all quaternions x. Note that the definition of the extension pair is unique up to 
units of 0, and that one has the relation 

nr(g Q ) = lcm(nr(g), nr(g )) = nr(g^) = \q a q a \ G N, 

which is crucial in the proof of the following Theorem. For the further characterisation of the 
CSLs of L, it is convenient to define the set 

(3) L(q) = {qx + xq \ x G 1} = <p + (ql) , 

which is a sublattice of L, compare Lemma [2. li Note that, due to I = I, one has L(q) = L{q). 

Theorem 4.4. Let q G I be admissible and primitive, and let q a = a q q be its extension. 
Then, the CSL defined by q is given by 

L n -T^qLq = L(q a ) , 

ml 

with L(q a ) defined as in Eq. ([3]). 

This explicit identification of the CSL provides access to the corresponding index 
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Theorem 4.5. If q £ I is an admissible primitive icosian, the rotation x i— > j^qxq is a 
coincidence isometry of L. Moreover, the corresponding coincidence index is 

E(q) = nr(q a ) = lcm(nr(g), nr(g)') , 

which is, with our above convention from Eq. fl2J), always an element of N. 

At this point, it is possible to determine the number of distinct coincidence rotations of 
index n, where n is a prime power. The rotations come in multiples of 120, the order of the 
rotation symmetry group of A4. 

Theorem 4.6. The number of coincidence rotations of a given index n is the same for the 
lattices A\,A± and L. Let 120/soc(^) be the number of coincidence rotations of index n. 
Then, fsoc{ n ) is a multiplicative arithmetic function, given at prime powers p r with r > 1 
by 

fsocipl = I jSp r -V +1 +P r - 1 - 2), p = ±1 (5) , 

p 2r + p 2r-2 ; p = ±2 (5) . 

Its Dirichlet series generating function reads 

' /socH _ Cjc(s-I) CWC(a-2) 



Dsoc{s) =J2 



n=l 



n s l + 5" s C(2s)C(2s-2) 



— l_ 5 2-s J_J_ (l-pl-s) (l-p2-s) J_J_ l- p 2- S 

p=±l(5) P=±2(5) 

_ ii 5 | 10 | 20 , 30 , 50 | 50 , 80 , 90 , 150 , 144 , 

2 s 3 s 4 s 5 s 6 s 7 s 8 s 9 s 10 s 11 s ' ' *' 

where £(s) is the Riemann zeta function and Ck( s ) denotes the Dedekind zeta function of the 
quadratic field K ; see Theorem 1 3. 81 for their explicit expressions. 

5. Outlook 

Due to the factorisations of icosians into irreducible elements, which are difficult to access, 
we have not yet found a way to approach the number of CSLs of A4 of index n systematically. 
Nevertheless, it is possible to derive the number of CSLs up to index 120, see [3]. We hope 
to report on some progress soon. 
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